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INTRODUCTION 

Environmental  conditions  at  the  time  of  ice  forma¬ 
tion  largely  determine  its  structure.  Ice  crystals  are 
uniaxial,  and  the  optic  axis  corresponds  to  the  c-axis. 
Random  c-axis  orientation  is  commonly  observed  near 
the  top  surface  of  newly  formed  sea  ice  in  the  Arctic. 
Once  acover  has  formed,  the  ice  structure  is  characterized 
by  long  vertical  columns  that  extend  downward  in  the 
growth  direction  of  the  ice  sheet,  a  result  of  quiescent, 
unidirectional  growth.  Under  these  conditions  a  selec¬ 
tive  growth  process  occurs,  and  the  c-axes  of  the  crys¬ 
tals  become  primarily  oriented  in  the  horizontal  plane  of 
the  ice  sheet  (Weeks  and  Ackley  1982).  In  the  presence 
of  a  predominant  current  direction,  strong  c-axis  align¬ 
ment  develops  in  the  direction  of  the  current  with 
generally  decreasing  scatter  as  ice  thickness  increases 
(Weeks  and  Gow  1978),  and  this  ice  structure  causes 
anisotropic  material  behavior  in  all  directions.  Uniaxial 
compression  data  on  first-year  sea  ice  (Richter-Menge 
et  al.  1987,  Wang  1979)  have  indicated  a  strong  depen¬ 
dence  of  peak  compressive  strength  on  c-axis  alignment 
and  on  the  angle  between  an  applied  load  and  the 
dominant  c-axis  direction.  Therefore,  to  interpret  data 
from  mechanical  property  tests  we  must  define  the 
relative  orientation  and  alignment  of  the  ice  fabric. 

The  techniques  used  in  the  analysis  of  ice  fabrics 
were  originally  developed  in  structural  petrology  (see, 
for  example,  Fairbaim  1 949,  Knopf  and  Ingerson  1 938, 
Turner  and  Weiss  1963).  Crystal  orientation  measure¬ 
ments  usually  involve  optical  measurements  of  the  c- 
axis  orientations.  Langway  ( 1 958)  describes  techniques 
for  obtaining  ice  crystal  c-axis  orientation  data  using  a 
Rigsby  universal  stage.  One  orientation  measurement  is 
made  for  each  ice  crystal  in  a  sample,  and  these  data  are 
plotted  on  a  Schmidt  equal  area  net  (Fig.  1 )  that  represents 
a  hemisphere  of  unit  radius.  The  points  on  the  net 
compose  an  orientation  diagram  that  depicts  the  relat  i  ve 
spatial  concentration  of  the  data,  and  this  representation 


reduces,  by  one,  the  dimension  of  the  hemisphere  and  of 
any  plane  or  line.  An  orientation  diagram  may  indicate 
a  random  or  a  patterned  structure,  depending  on  the 
dominant  features  of  the  diagram.  A  random  fabric 
refers  to  a  homogeneous  distribution  of  the  plotted 
points  and  represents  an  isotropic  material.  In  this 
configuration  there  is  an  equal  probability  of  finding 
points  in  equal  area  elements  anywhere  on  the  net.  In 
contrast,  the  most  significant  feature  of  an  anisotropic 
material  is  the  preferred  orientation  indicated  by  the 
grouping  of  points  on  the  net.  A  “girdle”  corresponds  to 
data  that  are  distributed  along  a  great  circle  of  the  net, 
indicating  a  preferred  planar  orientation  of  the  c-axes. 
The  pole  of  this  great  circle  is  termed  the  girdle  axis.  An 
area  of  highly  concentrated  data  points  (point  maxi¬ 
mum)  indicates  a  linear  preferred  orientation  of  the 
crystals  in  the  fabric.  The  statistical  significance  of  the 
orientation  diagram  increases  if  the  main  features 
(maxima/girdles)  are  reproducible  in  different  compa¬ 
rable  samples  from  the  same  homogeneous  body. 

Pearson  (1901)  used  statistical  arguments  to  develop 
the  equation  of  a  line  or  plane  that  provides  the  closest 
fit  to  points  in  space  when  all  variables  contain  error. 
The  solution  depends  on  knowledge  of  the  means, 
standard  deviations,  and  correlations  of  the  variables.  A 
significant  result  was  that  the  plane  of  best  fit  contains 
the  line  of  best  fit.  Watson  (1966)  presented  a  matrix  of 
sums  of  direction  cosines  of  vectors  representing  crystal 
orientation  in  a  Cartesian  coordinate  system.  The  rea¬ 
soning  presented  was  that  the  greatest  moment  of  inertia 
of  the  points  would  be  about  the  eigenvector  corre¬ 
sponding  to  the  minimum  eigenvalue  of  the  matrix. 
Mardia  (1972)  used  the  same  reasoning  to  obtain  this 
matrix  and  interpret  the  results.  Both  Watson  and  Mardia 
propose  distributions  for  the  data  on  the  sphere,  and 
develop  statistical  analyses  based  on  these  assump¬ 
tions.  Diggle  and  Fisher  (1985)  describe  a  program  that 
computes  these  eigenvalues  and  eigenvectors  and  quan¬ 
titatively  contours  spherical  data.  The  analysis  of  ice 


fabric  diagrams  has  been  largely  visual,  frequently 
based  on  approximate  data  concentration  contours  drawn 
on  the  net.  However,  Herron  and  Langway  ( 1 982)  ap¬ 
plied  the  eigenvalue/eigenvector  method  of  Mardia 
( 1972)  to  study  various  fabric  types,  including  small  cir¬ 
cle  girdles  and  multi-maxima  patterns.  The  results  were 
interpreted  qualitatively,  and  it  was  not  clear  whether  any 
fabrics  could  be  quantitatively  assessed  with  this  method. 

In  this  report  we  seek  a  plane  through  the  origin  that 
minimizes  the  sum  of  the  squared  normal  distances 
from  the  data,  and  obtain  the  dominant  c-axis  orienta¬ 
tion  in  this  plane.  Beginning  with  the  ice  orientation 
data,  a  detailed  derivation  from  simple  geometric  argu¬ 
ments  is  developed,  yielding  least-squares  equations 
that  minimize  the  orthogonal  distance  between  the  data 


and  a  best  line  and  plane.  The  eigenvalue  problem  that 
results  is  the  same  as  that  obtained  by  Pearson,  Watson 
and  Mardia.  We  identify  an  implicit  assumption  in  this 
method  of  equal  measurement  uncertainty  in  each  coor¬ 
dinate  at  all  points  on  the  sphere.  Normalized  eigenval¬ 
ues  provide  quantitative  measures  of  physical  distance 
of  the  data  from  the  plane  and  line,  specifying  the 
directional  characteristics  of  the  c-axes  of  crystals  in  a 
sample.  Mean  angular  measures  of  variability  are  also 
developed.  The  development  clearly  indicates  the  fab¬ 
rics  that  are  well-described  by  the  eigenvalue/eigenvec¬ 
tor  method,  and  provides  a  framework  for  developing 
related  methods  that  may  be  useful  for  quantifying  other 
fabric  types.  We  demonstrate  the  capabilities  of  the 
analysis  on  data  sets  representing  several  samples  of 
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Figure  1 .  Schmidt  equal  area  net. 


first-year  sea  ice.  The  data  and  the  results  of  the  analysis 
are  viewed  on  special  Schmidt  nets  that  represent  data 
hemispheres  defined  by  the  best  plane  and  the  predomi¬ 
nant  basal  plane,  in  addition  to  the  Ay-plane. 

BACKGROUND 

Measurement  of  the  c-axis  orientation  of  an  ice 
crystal  with  a  Rigsby  universal  stage  provides  an  azi¬ 
muth  angle,  an  inclination  angle,  the  direction  of  in¬ 
clination  as  right  or  left,  and  the  type  of  measurement  as 
polar  or  equatorial.  From  these  data  we  will  obtain  apair 
of  angles  that  define  the  orientation.  The  azimuth  angle 
measures  the  rotation  of  the  crystal  about  the  z-axis  that 
places  the  c-axis  in  the  xz-plane.  The  zero  azimuth  can 
be  chosen  arbitrarily,  depending  on  the  requirements  of 
the  analysis.  With  the  c-axis  in  thexz-plane,  the  inclina¬ 
tion  angle  is  the  angle  of  tilt  with  the  v-axis.  An  inclina¬ 
tion  angle  of  zero  is  defined  as  the  angle  where  the  c-axis 
of  the  crystal  is  parallel  to  the  x-axis.  It  is  found  by 
rotating  the  crystal  about  the  y-axis  either  to  the  right 
(clockwise)  or  the  left  (counterclockwise).  The  optic 
axis  of  a  crystal  can  be  aligned  parallel  to  the  z-axis, 
termed  a  polar  measurement,  orparallel  to  the  x-axis,  an 
equatorial  measurement.  The  inclination  angle  measured 
with  the  Rigsby  stage  must  be  corrected  for  the  optical 
error  caused  by  the  difference  in  the  refractive  indices 
of  air  and  ice.  The  relationship  between  the  corrected  in¬ 
clination,  /,  and  the  measured  inclination,  /m,  depends 
on  the  measurement  type  and  can  be  expressed  (Kamb 
1962)  as 


circle  planes  passing  through  the  y-axis  with  the  unit 
sphere.  These  arcs  are  meridians  of  the  net  that  corre¬ 
spond  to  inclination  angles  between  -90°  and  +90°, 
representing  the  angle  between  the  lower  half  of  the  yz- 
plane  and  the  plane  of  the  arc.  The  labels  on  the  x-axis 
correspond  to  polar  or  equatorial  crystal  orientation 
with  right  or  left  tilt.  The  set  of  arcs  that  cross  the  y-axis 
are  parallels  of  the  net,  and  represent  the  intersection  of 
planes  parallel  to  thexz-plane  with  the  sphere.  The  points 
of  intersection  of  these  arcs  with  the  perimeter  of  the  net 
represent  the  azimuth  angle,  and  the  notation  on  the  net 
is  in  degrees  of  angle.  Additional  details  that  are  needed 
for  plotting  crystal  orientation  data  on  the  Schmidt  net 
are  given  by  Langway  (1958). 

DATA  TRANSFORMATION  TO  CARTESIAN 
COORDINATES 

Our  analysis  represents  the  c-axis  of  each  ice  crystal 
as  a  unit  vector  from  the  origin  of  a  three-dimensional 
coordinate  system  in  the  half-space  below  thexy-plane, 
yielding  an  array  of  points  on  the  surface  of  a  hemi¬ 
sphere  of  unit  radius.  In  this  section  we  obtain  three- 
dimensional  Cartesian  coordinates  on  the  unit  hemi¬ 
sphere  (Fig.  2)  for  each  crystal.  The  initial  step  in 
finding  these  Cartesian  coordinates  is  to  represent  each 
point  in  spherical  coordinates  (p,0,<j>).  Because  each 
crystal  is  represented  by  a  unit  vector,  the  radius  p  =  1 
for  all  the  data.  The  angle  0,  measured  from  the  positive 
x-axis,  is  positive  in  the  counterclockwise  direction,  and 
can  be  obtained  from  the  measured  azimuth  Az  as 


1.04 /m 

;  equatorial 

sin"1  I — ! — sin  /m) 

,  11.31  / 

;  polar 

Taken  together,  the  azimuth  and  inclination  mea¬ 
surements  define  the  orientation  of  the  optic  axis  of  an 
ice  ciystal  in  three-dimensional  space.  The  line  repre¬ 
senting  each  crystal  in  the  sample  is  plotted  through  the 
origin  of  a  sphere  of  unit  radius.  Each  of  these  lines 
intersects  the  surface  of  the  sphere  at  one  point  in  each 
hemisphere.  The  surface  of  the  lower  hemisphere  is 
traditionally  represented  in  two  dimensions  with  the 
Schmidt  net  (Fig.  1 ).  The  Schmidt  net  is  also  known  as 
the  equal-area  net  because  a  unit  area  in  any  position  on 
the  net  corresponds  to  a  unit  area  on  the  hemisphere.  The 
crystal  orientation  data  plotted  on  this  net  depict  the 
relative  spatial  concentration  of  the  measured  data.  The 
Schmidt  net  is  composed  of  two  sets  of  arcs.  The  set  that 
crosses  the  x-axis  represents  the  intersection  of  great 


a  Az  ;  polar-right,  equatorial-left  1 

0=  (2) 

[Az  +  180°  ;  polar-left,  equatorial-right  ) 

where  Az  is  in  degrees.  Measurements  made  from  the 
negative  x-axis  are  adjusted  by  1 80°  in  eq  2. 

The  angle  <j)  is  measured  from  the  positive  z-axis.  The 
inclination  angle  was  measured  from  the  negative  z-axis 
for  polar  crystals  and  from  the  xy-plane  for  equatorial 
crystals,  and  <|)  is  determined  as 


1 1 80°  -  /  :  polar  j 
\  90°  +  /  ;  equatorial  I 


(3) 


With  the  spherical  coordinates  specified,  the  equiva¬ 
lent  Cartesian  coordinates  are  obtained  directly  as 

x  =  p  sin<j>  cos0 

y  =  p  sin<()  sin9  (4) 

z  =  p  cos<|> 
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z 
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Figure  2.  Sketch  of  unit  vectors  1 ,2,..., m,..fJ  representing  the  c-axis  orien¬ 
tations  of  ice  crystals  in  a  sample.  The  hemisphere,  z  <  0,  is  shown  by  con¬ 
vention. 


MEASUREMENT  ERRORS  AND 
DATA  FITTING 

The  relative  uncertainties  in  x,  y,  and  z  caused  by 
measurement  error  must  be  understood  in  order  to 
identify  an  appropriate  method  for  fining  the  ice  fabric 
data.  Langway  (1958)  lists  several  sources  of  error  in 
the  measurement  of  c-axis  orientation,  and  indicates 
that  the  errors  in  aximuth  and  inclination  should  each  be 
less  than  5°.  Taking  the  errors  in  Az  and  /  as  random  and 
of  comparable  magnitude,  we  observe  from  eq  2  and  3 
that  the  same  statement  holds  for  errors  in  0  and  <|>.  For 
an  error  of  P  in  these  angles  we  have  from  eq  4  that 

x  =  sin(<t>  ±  P)  cos(0  ±  P) 

y  =  sin(<p  ±  P)  sin(0  ±  P)  (5) 

z  =  cos(<|»  ±  P) 

Expanding  eq  5  and  grouping  terms  according  to  the 
error  we  obtain 

x  =  sirup  cos0  cos2P  +  sinp  (A  cosP  +  B  sinP) 
y  =  sinp  sin0  cos2P  +  sinp  (C  cosP  +  D  sinP)  (6) 
z  =  cosp  cosP  +sinp  (E) 

where  A,B,C,D,E  are  composed  of  sines  and  cosines 
of  p  and  0.  The  relationships  given  in  eq  4  are  contained 
in  eq  6,  but  they  have  been  modified  by  errors  of  p.  The 


effect  of  measurement  errors  on  the  uncertainty  in  x,  y, 
and  z  is  variable.  An  error  in  0  only  affects  x  and  y.  The 
resulting  uncertainty  in  each  coordinate  is  periodic, 
depending  on  0,  and  out  of  phase  with  the  other.  The 
amplitude  of  these  errors  approaches  zero  near  the  pole 
and  a  maximum  at  the  equator.  Errors  in  p  affect  all  three 
Cartesian  coordinates.  The  uncertainty  in  z  is  larger  near 
the  equator  and  smaller  near  the  pole  than  the  larger  of 
x  and y.  Uncertainties  in  x  and  y  resulting  from  errors  in 
p  again  vary  individually  with  0,  displaying  maximum 
ampl  itudes  near  the  pole  and  approaching  zero  near  the 
equator. 

Classical  least-squares  methods  require  a  dependent 
variable,  but  the  coordinates  of  each  c-axis  are  all 
independent.  Reed  (1989)  presented  a  method  for  fit¬ 
ting  a  line  to  points  in  the  plane  when  both  coordinates 
of  these  points  are  independent  and  uncertain  due  to 
measurement  error.  The  method  allows  the  errors  in  the 
fit  to  be  weighted  according  to  the  relative  uncertainties 
in  the  measurement  of  x  andy.  The  measurement  errors 
discussed  above  could  be  considered  by  extending  this 
method  to  three  dimensions.  However,  we  note  for 
small  P  that  cosP = cos2p  =  1  and  sinp = 0,  indicating  that 
errors  in  x,  y,  and  z  resulting  from  measurement  errors 
are  relatively  small.  As  a  first  approximation  we  will 
choose  equal  weighting  in  each  direction  at  every  point, 
and  the  best  fit  line  and  plane  will  minimize  the  perpen¬ 
dicular  distances  from  the  data. 
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DETERMINATION  OF  THE  BEST  PLANE 
BY  ORTHOGONAL  LEAST-SQUARES 

The  unit  vectors  representing  the  c-axis  orientation 
of  each  crystal  in  a  sample  have  a  common  point  at  the 
origin  of  the  unit  sphere.  The  problem  we  consider  in 
this  section  is  to  find  the  plane  of  best  fit  to  fabric  data 
that  contains  the  origin,  and  to  provide  quantitative 
measures  of  the  quality  of  the  fit.  The  form  of  the 
equation  of  a  plane  through  the  origin  is 

f(x,y,z)  =  Ax  + By +  Cz-0  (7) 

We  choose  the  function  F(x,y,z;c  ( ,c2,c3)  with  this  same 
form  and  depending  linearly  on  parameters  c , ,  c2,  as 

3 

F(x,y,z\c\,  C2,  C3)  =  ci<J>!  +  C2<J>2  +  c3<J>3  =  £  Cjfo  =  0 

i  =  1 

(8) 

where  =x,  <p2  =y,  <t»3 = 2  are  a  specified  set  of  mutually 
orthogonal  functions  and  the  Cj  are  unknowns  to  be 
determined.  The  pole  of  a  plane  is  the  point  P  of  inter¬ 
section  with  the  hemisphere  of  a  line  normal  to  the  plane 
through  the  origin.  The  unit  normal  n  to  the  plane  of  best 
fit  to  the  data  has  the  form 

n  =ci  ii  +  c2i2  + c3i3  =  Cjij  (9) 

where  i,,  i2  and  i3  are  unit  vectors  in  the  x,  y  and  z 
directions,  respectively,  and  repeated  indices  indicate 
summation. 

A  unit  vector  in  3-dimensional  space  represents  each 
crystal  in  a  sample,  and  the  total  number  of  crystals  N 


used  to  fit  the  plane  should  be  much  larger  than  3.  For 
a  data  point  m  located  at  (-vm,ym,rm)  we  can  evaluate 
C1-C2’C3>  =  fm  with  eq  8  as 

^ m  ~  c  I  -'in  +  c3“m  (10) 

The  unit  vector  representing  the  c-axis  of  the  mth  crys¬ 
tal  (Fig.  3)  is 

um  =  xm  i  1  +  ym  *  2  +  Zm*  3  =  ^  1  1  + 

<t>2m*2  +  4>3m*3  =  ^jm'j  (I*) 

The  vector  wm  is  the  projection  of  um  onto  n,  rep¬ 
resenting  the  normal  vector  from  the  plane  to  the  point 
(■*nvym’zm) 

wm  =  (um  -n)n  (12) 

where  (um  •  n)  is  a  scalar  product  between  unit  vectors. 
Throughout  this  development  the  notation  (~)  over  a 
vector  indicates  that  it  does  not  have  unit  length.  The 
normal  distance  dm  from  point  m  to  the  plane  is 

4m  =  (um  •  n)  =  C,Jfm  +  QJm  +  C32m  ( ]  3) 

and  we  observe  that 

KJ  =  IWml  S  1  (14) 

wm  =  4nn 

The  sign  of  dm  distinguishes  distances  on  opposite  sides 
of  the  plane. 

The  orthogonal  projectionof  um  in  the  best  fit  plane 
is  represented  by  the  vector  vm,  and 


Figure  3 .  Sketch  of  the  unit  vector  um  representing  the  mth  crystal  in  an  ice  sample  and 
its  projections  vm  onto  the  best-fit  plane  or  the  basal  plane,  and  wm  onto  norc  the 
unit  normal  vector  to  the  corresponding  plane.  The  angles  between  um  and  the  best 
(plane,  line)  are  (<xm,  &m),  respectively. 
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vm=um-  wm  =  xmi  |  +  ymi2  +  zmi3  =  (15) 

This  vectorjoins  the  origin  and  the  point  a  m,  vm,  rm,the 
orthogonal  projection  of  (.vm  ,ym  ,rm )  in  the  plane.  The 
vector  vm  is  unique  and  the  best  approximation  of  um  in 
the  plane.  An  equivalent  statement  is  that  the  closest 
point  in  the  plane  to  (VV’W  is  Am.  ym.  :m.  From  eq 
15  we  can  determine  .vm,  vm,  :m  as 

£m  =  'm  ' 

vm=ym-4nf2  for  m  =  1,2 N  (16) 

-in  =  -m  ’  4nO 

or  ~ 

0im  =  0im  -  <Wi  for/  =1,2.3  (17) 

where  <t>im  and  0^  represent  elements  of  iV-dimensional 
vectors  <t>|  and  (t»j,  respectively. 

We  want  tochoose  the  c-  values  that  speci  fy  the  plane 
through  the  origin  with  normal  distances  Wml  between 
the  points  representing  the N crystals  in  a  sample  and  the 
plane  that  are  as  small  as  possible.  The  A/-vectorr/  =  (<£, 
d2, ....  dN)T  represents  the  individual  normal  distances 
from  the  data  to  the  plane,  where  T  indicates  the  trans¬ 
pose.  Minimum  normal  distance  to  the  plane  is  equiva¬ 
lent  to  maximum  length  of  v  m,  the  projection  of  um  onto 
the  plane.  We  will  seek  a  least-squares  fit  and  define  £ 
as  a  function  of  the  unknown  coefficients 

N  N  3  _  , 

£(q.Q.Q)=  X  Ivmp=  I  (18) 

m  =  1  m  =  I  /  =  I 

Because  we  consider  data  projected  onto  a  plane,  the 
hemisphere  in  which  the  data  appears  is  arbitrary.  The 
sum  of  the  squares  of  the  lengths  of  vm  will  be  maxi¬ 
mized  and  the  sum  of  the  squared  distances  from  the 
plane  dm  will  be  minimized  where  the  gradient  of  £ 
vanishes 

VE(c].c2,c*i)  =  0  (|9) 

with  Cj*  indicating  coefficients  of  the  plane  of  best  fit. 

We  now  perform  the  differentiation  indicated  in  eq 
1 9  with  respect  to  Cj,  where  each  choice  of /yields  a  scalar 
equation  of  the  form 


where  5(J  is  the  Kronecker delta.  Vi  » observe  from  eq  1 3 
that  dm  =  ek  <j>km,  and  differentiation  of  this  equation 
yields 

~  0  km  Sjk  =  0jm  (22) 

dfj 

Inserting  eq  17,  21  and  22  into  eq  20  we  obtain  the 
normal  equations  as 


or  expanding  and  rearranging  as 

0jn/4n  ( 1  “  )  +  0im  0jm  f  i  =  ^ m  c’j  (23) 

Finally,  because  n  is  a  unit  vector  and  /  is  a  summation 
index,  (1  -  c7  )  =  0  and  eq  23  becomes 

<t>im0jm<i  =  (j  (24> 

In  eq  24  we  observe  that  dm~  =  X,  a  constant,  yielding 
an  eigenvalue  problem, 

Ac  =  kc  (25) 

where  each  term  a,j  of  the  3  x  3  symmetric  matrix  A  is 

oij  =  0im0jm  =$1-  0j  (26) 

The  eigenvalue  k  is  the  sum  of  the  squared  normal 
distances  of  the  data  from  the  corresponding  plane. 
Equation  25  in  homogeneous  form. 

(A  -  A/)  <  =  0 .  (27) 

indicates  that  nontrivial  solutions  exist  if  and  only  if 

dct{A  -  kf)  =  0  (28) 

The  determinant  given  in  eq  28  yields  a  cubic  equation 
called  the  characteristic  polynomial 

A2  +  pk~  +  qk+  /=  0  (29) 

Following  Beyer  (1987)  we  obtain  the  solution  as 
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m  =  2  V  -aji 

a-  (1/3)  (3r/  —  p 2) 

b  =  (l/27)(2p3-  9pq+  27 r) 


As  p  and  r  are  real,  the  eigenvalues  will  be  real  and 
tv  met  if 


(31) 


The  spectral  theorem  (Shields  1968)  states  that  be¬ 
cause  A  is  a  real  symmetric  matrix  it  is  similar  to  a 
diagonal  matrix  B  composed  of  the  eigenvalues  of/f,  and 
therefore  the  eigenvalues  of  A  are  real.  Similar  matrices 
have  the  same  trace  and  the  same  determinant.  The  trace 
of  a  matrix  is  the  sum  of  the  elements  on  the  principal 
diagonal.  The  elements  o,j,  a22  and  of  matrix  A 
represent  the  sum  of  the  squares  of  the  distances  be¬ 
tween  the  data  and  the  three  planes  defined  by  the 
coordinate  axes.  As  each  crystal  is  represented  by  a  unit 
vector, 

3  3  3 

tr(A)  =  X  °ii  =N  =  ti\B)  =  X  &ii  =  X  *■'  (32) 

i  =  I  i= 1  i =  1 

indicating  that  the  sum  of  the  eigenvalues  is  N,  the  total 
number  of  crystals  in  the  sample.  The  determinants  of  A 
and  B  are  the  product  of  the  eigenvalues.  The  eigen- 
valuesofA  are  nonnegative  iM  is  positive  semidefinite, 
that  is.  .vtA.v  >  0  for  any  x.  Large  diagonal  elements 
relative  to  those  off  the  main  diagonal  are  characteristic 
of  a  positive  semidefinite  matrix.  Since  the  eigenvalues 
represent  sums  of  squared  distances,  X(  >  0  and  matrix 
A  is  positive  semi  definite.  We  will  designate  the  eigen¬ 
values  in  increasing  order  according  to  magnitude  as  X  ( 
^  Xo  <  X2,  and  define  normalized  eigenvalues  as 


(33) 


The  normalized  eigenvalues  give  the  mean  squared 
normal  distance  between  the  points  on  the  unit  sphere 
defined  by  the  c-axis  vector  of  the  individual  crystals 
and  the  plane  normal  to  the  corresponding  eigenvector. 
These  eigenvalues  provide  a  measure  of  the  fit  that  is 
equivalent  to  the  variance  in  classical  dependent  vari¬ 
able  least-squares  methods. 

Eigenvectors  of  a  real  symmetric  matrix  correspond¬ 
ing  todifferent  eigenvalues  are  orthogonal ,  and  because 
the  eigenvalues  are  real,  the  eigenvectors  can  be  taken 
to  be  real .  The  vector  v  is  an  eigenvector  for. 4  belonging 
to  the  eigenvalue  X  if 


Ay  =  Xv  and  v  *  0 


(34) 


basis  in  three-dimensional  space.  Each  eigenvector 
represents  the  unit  normal  to  a  plane,  and  the  corre¬ 
sponding  eigenvalue  gives  the  sum  of  the  squared 
normal  deviations  of  the  data  from  that  plane.  The 
minimum  eigenvalue  defines  the  plane  of  best  least- 
squares  fit  to  the  data,  and  the  higher  eigenvalues  are 
associated  with  the  remaining  mutually  orthogonal 
planes  through  the  origin.  With  the  origin  fixed,  the 
eigenvector  basis  represents  a  coordinate  system  that  is 
rotated  relative  to  the  coordinate  axes. 

The  eigenvectors  written  in  columns  form  the  matrix 
P.  The  elements  of  P  are  the  direction  cosines  between 
each  eigenvector  and  the  coordinate  axes.  The  angle  a- 
between  the  eigenvector  Vj  and  the  axis  ij  is 

cty  =  cos'1  (Vj  •  ij)  =  cos-1  (p|j)  (35) 

Because  the  columns  of  P  are  orthonormal.  P  is  an  or¬ 
thogonal  matrix,  P~ 1  =/>T  is  also  orthogonal,  and  Jet  P 
=  +  1 .  Matrices  A  and  B  are  related  through  P  as 

A  =  PBPJ  (36) 

representing  a  singular  value  decomposition  of  A.  The 
diagonal  elements  of  B  are  the  singular  values  Oj  as  well 
as  the  eigenvalues  of  A .  With  singular  values  ordered  by 
their  magnitude  in  the  same  way  as  the  eigenvalues,  the 
condition  number  of  matrix  A  of  full  rank  is 

cond  (A)  =5i  =  (37) 

Oi  X, 

Condition  numbers  »1  indicate  that  A  is  nearly  sin¬ 
gular. 

The  sum  of  the  squared  normal  distances  between 
the  data  and  the  best  plane,  given  by  Xj,  provides  a 
measure  of  the  planar  structure  of  an  ice  sample.  Values 
of  X|  or  Xj  approaching  zero  indicate  an  increas¬ 
ingly  planar  ice  structure.  If  X(  =  0.  then  cond(  A)  =  °° 
and  the  data  are  perfectly  planar.  A  visual  representa¬ 
tion  of  the  error  is  obtained  from  the  angle  am  between 
an  individual  crystal  and  its  projection  in  the  best  plane 
(Fig.  3). 

am  =  |7t/2-  cos"1  (um  •  n)  | 

(38) 

The  absolute  value  in  eq  38  is  needed  if  the  angle 
between  um  and  n  is  greater  than  rc/2.  The  average  an¬ 
gular  deviation  a  between  the  data  and  the  best  plane  is 
a  parameter  we  term  the  planar  spread  that  can  be 
readily  detenmined  and  understood: 


a 


m 


The  lengths  of  these  eigenvectors  are  arbitrary,  and  we 
normalize  them  to  unit  length  to  obtain  an  orthonormal 


(39) 


1 


A  small  planar  spread  indicates  a  small  mean  angle 
between  the  crystals  in  a  sample  and  the  plane. 


n  = 


(if + + 1 ) 


1/2 


'I 


1  i  +  ('-> 


•2-  '3) 


t45) 


DEVELOPMENT  OF  DEPENDENT  The  normal  distances  between  the  data  and  the  plane 

VARIABLE  LEAST-SQUARES  SOLUTIONS  represented  by  n  can  then  be  determined  using  eq  13  as 


We  will  now  parallel  the  development  of  the  previ¬ 
ous  section  to  obtain  a  set  of  classical  dependent  vari¬ 
able  least-squares  solutions  to  the  problem  of  fitting  a 
plane  to  three-dimensional  data.  This  standard  method 
seeks  a  minimum  error  in  the  dependent  variable  and  is 
more  commonly  available  than  the  orthogonal  method, 
but  the  results  are  sensitive  to  the  choice  of  dependent 
variable.  Arbitrarily  selecting  z  as  the  dependent  vari¬ 
able,  we  rewrite  eq  8  as 


:  =  EfA-.y.t  =  c,  <(),  +  c2  <|>2  (40) 

where  4> ,  =  x  and  <p-,  =  y.  The  coefficients  c,  and  c~,  are 
chosen  so  that  the  deviations  dm  are  as  small  as  possible. 


^m  =  rm-Hvm.ym:fV‘'2>'  (41> 

in  which  rm  is  zm,  the  r-value  of  the  mth  point,  and 
F(.xm,ym:c  t  ,c2)  is  the  corresponding  r-value  in  the  plane 
of  best  fit.  Again,  the  vector  d  contains  all  the  deviations 
dm.  and  we  seek  a  least-squares  fit  to  the  data  with 


N  2 

£ (n.  q)  =  X  [rm-£ (vm.  :ym;  q,  q)] 

m  =1 


(42) 


Taking  the  gradient  of  E  and  setting  it  to  zero  yields  the 
normal  equations 


[fin  £  (-'m>  3W  Q)]  —  0  *  i  —  \, 2 

m  =  I 

(43) 

or 

dT  •  =  0 


lwnJ=(u,n>n)  JC  |  Vnl  +  (lV?-  ;m)- 
(<-?  +  r?  +  l),/2 


(46) 


for  m  =  1 ,2 . N 


The  sum  of  the  squared  normal  distances  given  in  eq  46 
is  the  same  measure  as  the  minimum  eigenvalue  ob¬ 
tained  in  the  previous  section.  The  planes  that  minimize 
the  squared  .v-deviations  and  squared  y-deviations  are 
obtained  by  interchanging  the  roles  of  <J>,,  <j>2  and  r 
above,  and  repeating  the  analysis. 


ALIGNMENT  OF  THE  C-AXES 

An  area  of  highly  concentrated  c-axis  data  on  the 
Schmidt  net  indicates  a  linear  preferred  orientation  of 
the  ice  crystals  of  a  sample,  and  suggests  the  need  to 
determine  the  predominant  optic  axis  orientation  of  the 
fabric.  We  will  locate  this  linear  orientation  by  follow¬ 
ing  a  development  parallel  to  that  used  to  determine  the 
best  plane.  The  unit  vector  c  represents  the  unknown 
preferred  c-axis  orientation  of  the  crystals  in  a  sample, 
and  is  expressed  as 


c  =  ej  i  1  +  c,  i  2  +  Cj  i  3 

where  the  primes  distinguish  these  coefficients  from 
those  of  the  unit  vector  n  given  in  eq  10.  The  plane 
through  the  origin  that  is  normal  to  c  represents  the 
predominant  basal  plane  orientation  of  the  crystals  in 
the  sample  and  is  described  by 


indicating  that  the  vector  of  deviations  is  normal  to  each 
vector  ty.  Now,  inserting  eq  40  and  rewriting  the  normal 
equations  in  vector  form  we  obtain  a  pair  of  linear 
equations: 

[<t>I  *  <t>i]  =  rT  *  <j>i  ‘  =  L2  (44) 

The  vector  normal  to  the  plane  that  minimizes  the 
sum  of  the  squared  z-distances  with  the  data  is 

"  =  c,  i,+  c2  i2  -  h 
or  in  the  form  of  a  unit  vector. 


F  (  v’-v -  ri  ’ c2’ c3^  =  c\x  s  c2>’  +  c3:  =  0  (48) 

As  before,  the  unit  vector  um  represents  the  mth  crystal 
orientation  and  intersects  the  unit  hemisphere  at 
(.Ym,ym,zm).  In  order  to  take  advantage  of  the  detail 
given  in  the  previous  development,  the  vector  projec¬ 
tions  of  Up,  onto  the  plane  and  its  normal  vector  are  again 
vm  and  wm,  respectively.  Then,  with  c  replacing  n  and 
a  different  normal  plane.  Figure  3  represents  our  present 
condition. 

When  searching  for  the  best  plane  we  sought  to 
minimize  the  squared  normal  distance  lwml.  However, 
the  closest  representation  of  um  by  the  unit  vector  c 
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Figure  4.  Unit  vector  c  representing  the  linear  preferred  c-a.xis  orientation  of  the 
ice  fabric,  and  the  long  axis  of  the  columnar  crystals  of  an  ice  sample  are 
represented  by  n.  The  unit  vectors  n  and  s  provide  the  predominant  orientation 
of  the  basal  plane  of  the  sample.  A  unit  vector  a  represents  the  direction  of  toad 
application  on  the  sample,  a  is  the  projection  of  this  vector  onto  the  basal  plane, 
ac  is  the  angle  between  the  applied  load  and  c,  and  a,  is  the  angle  between  the 
load  and  n. 


requires  that  we  maximize  the  sum  of  the  squares  of  the 
lengths  of  w  m.  The  predominant  basal  plane  orientation 
is  then  the  plane  of  maximum  squared  normal  error  with 
the  data.  Each  step  in  the  previous  development  applies 
except  that  c;  is  replaced  by  d.  Note  that  the  matrix  A 
given  in  eq  26  is  unchanged  because  it  depends  only  on 
the  coefficients  of  the  intersections  of  the  individual 
unit  vectors  um  with  the  unit  hemisphere.  Therefore,  we 
are  solving  the  same  eigenvalue  problem  as  before.  The 
eigenvectors  obtained  in  eq  34  are  orthogonal,  and  c  is 
contained  in  the  best  plane.  Kass  (1989)  has  shown  for 
the  general  case  that  spaces  of  closest  fit  are  nested. 
Every  p-dimensional  subspace  of  closest  fit  lies  in  one 
of  dimension  p  + 1 .  The  third  eigenvector  s  together  with 
c  form  the  plane  of  best  fit  to  the  c-axes,  ands  and  n  form 
the  predominant  basal  plane  (Fig.  4).  The  direction 
cosines  of  the  angles  between  the  preferred  1  inear  orien¬ 
tation  and  the  coordinate  axes  are  given  in  matrix  P,  and 
the  angles  can  be  obtained  from  eq  35. 

The  mean  squared  normal  distance  between  the  data 
andthevectorcis  X',+  X'2orl  -  Xy  Again  using  angles 
to  visualize  error,  the  angle  8m  between  an  individual 
(crystal)  vector  and  its  projection  onto  c  (Fig.  3)  is  ob¬ 
tained  as 

8m  =  cos-1  (um  *c)  (49) 


Linear  spread  between  these  vectors  is  al  ways  <  90°. 
If  8m  obtained  with  eq  49  is  >  9(F,  8m  is  replaced  by  its 
supplement.  The  average  angle  8  between  the  data  and 
c  is 

*  =  jj  I  8"  (50) 

The  parameter  8  is  termed  the  linear  spread,  and  a  small 
8  indicates  that  the  angles  between  the  preferred  orien¬ 
tation  and  the  data  are  also  small. 

As  a  group,  the  eigenvalues  provide  measures  of  the 
structure  of  the  ice  sample.  An  eigenvalue  of  zero 
occurs  if  the  data  are  2-dimensional,  and  a  pair  of  zero 
eigenvalues  represent  1  -dimensional  data.  More  gener¬ 
ally,  small  X'f  «  X'2  together  with  large  X'3  occurs  with 
aligned  data,  and  small  X.',  with  significantly  larger  X'2 
and  X3  indicates  planar  data.  If  in  addition  to  small  X')t 
X'2  =  X'3  the  crystal  orientations  distribute  symmetri¬ 
cally  about  a  great  circle  of  the  sphere.  With  X',  =  X'2 
=  X'3  every  set  of  orthonormal  vectors  will  serve  as  the 
eigenvectors.  The  data  are  maximally  dispersed  with 
respect  to  lines  and  planes,  and  the  material  is  isotropic. 

If  the  c-axes  in  a  given  sample  are  sufficiently 
aligned,  the  mechanical  properties  of  the  ice  will  be 
affected.  If  a  is  a  unit  vector  in  the  direction  of  an 
applied  force,  the  angle  ctc  between  the  force  and  the 
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dominant  c-axis  direction  is 

oc  =  cos-1  (o  •  c)  (51) 

and  the  complement  of  cc  is  the  angle  between  the  load 
and  the  basal  plane  (Fig.  4).  For  a  columnar  ice  sample 
n  gives  the  predominant  direction  of  crystal  elongation 
and  growth.  The  angle  az  between  the  load  and  the 
vector  n  is  obtained  as 

az=cos~l  (o»  n)  (52) 


SCHMIDT  NET  REPRESENTATIONS 

The  results  of  the  analysis  of  an  ice  sample  are 
presented  together  with  the  data  on  the  Schmidt  net .  The 
initial  step  in  finding  the  Cartesian  net  coordinates  is  to 
obtain  polar  coordinates  (r,0)  on  the  net  for  each  point. 
The  polar  angle  0  on  the  net  is  the  same  as  the  spherical 
angle  0  given  in  eq  2.  The  spherical  angle  <j>  given  in  eq 
3  is  used  to  find  the  radial  coordinate.  The  supplement 
of  <J>,  measured  from  the  negative  z-axis,  has  the  same 
sine  as  <j>,  and  the  distance  from  the  origin  of  the  Schmidt 
net  r  is  (Knopf  and  Ingerson  1938) 

r  =  (2a)sin(<)>/2)  (53) 

where  a  =  0.7071 1/?,  and  R  is  the  radius  of  the  net.  The 
equivalent  Cartesian  coordinates  are  then  obtained  di¬ 
rectly  as 

x  =  r  cos0 

(54) 

y  =  r  sin0 

Planes  that  pass  through  the  origin  intersect  the  unit 
sphere  as  a  great  circle.  The  part  of  the  great  circle  of  the 
best  fit  plane  below  the  Ay-plane  is  drawn  on  the  Schmidt 
net.  With  p  =  1 ,  eq  4  is  substituted  into  eq  8  and  solved 
for  <|>  to  obtain 

<j>  =  tan-1  ( - =£2 - -j  (55) 

\ci  cos0  +  02  sin 0/ 

where  0  is  incremented  in  arbitrary  steps  from  OP  to  360°. 
The  requirement  of  z  £  0  in  eq  4  identifies  the  points  of 
intersection  in  the  lower  hemisphere.  We  obtain  polar 
coordinates  on  the  net  from  these  spherical  coordinates, 
and  then  eq  54  yields  a  discrete  representation  in  Carte¬ 
sian  coordinates  of  the  great  circle  on  the  net.  Experi¬ 
ence  indicates  that  0.5°  increments  of  0  yield  a  smooth 
curve  on  the  net. 

The  unit  normal  n  to  the  plane  of  best  fit  from  the 


origin  is  the  pole  vector.  It  intersects  the  unit  hemi¬ 
sphere  at  the  pole  point  P  with  coordinates  (.vp,>'p,rp)  or 
(C|  .Cj.cj).  Inverting  eq  4  we  obtain  the  spherical  coor¬ 
dinates  of  this  point  as 

<(>  =  co^1  (c3)  (56) 

9=,an-'(a) 

There  are  two  values  of  0  in  the  range  0°  <  0  <  360°  that 
have  the  same  tangent.  If  the  value  of  a-  =  Cj  is  negative, 
0  falls  in  the  2nd  or  3rd  quadrant  and  the  calculated 
angle  is  adjusted  by  adding  1 80°.  Again,  the  polar  co¬ 
ordinate  r  on  the  net  is  obtained  from  eq  53,  and  the 
Cartesian  coordinates  of  the  pole  on  the  net  are  found 
from  eq  54.  The  determination  of  the  Schmidt  net 
coordinates  of  c,  representing  the  intersection  of  the 
preferred  linear  orientation  and  the  unit  hemisphere, 
follows  the  same  procedure  as  for  the  pole  vector  with 
(c1(c2,c3)  in  eq  56  replaced  by  (c.  ,c'2,c3).  The  intercept 
of  the  linear  orientation  vector  must  fall  on  that  of  the 
great  circle  of  best  planar  fit. 

During  field  ice  coring  and  thin  section  preparation, 
a  sample  intended  as  horizontal  may  deviate  by  an  angle 
of  several  degrees.  Sea  ice  with  horizontal  c-axis  align¬ 
ment  will  then  be  represented  by  a  best  plane  at  this 
angle.  The  linear  dimensions  on  the  periphery  of  the 
Schmidt  net  are  distorted,  and  it  is  difficult  to  judge 
normal  distances  between  the  data  in  this  region  and  the 
trace  of  the  best  plane.  Points  that  appear  near  the 
perimeter  and  directly  across  the  net  from  each  other 
represent  crystals  with  close  planar  alignment.  For 
these  reasons  it  is  frequently  advantageous  to  view  the 
data  on  Schmidt  nets  drawn  on  alternative  planes. 
Viewing  the  data  on  a  Schmidt  net  drawn  in  the  best  fit 
plane  with  n  vertical  eliminates  sample  preparation  er¬ 
ror  for  horizontally  aligned  sea  ice.  The  great  circle  of 
the  best  plane  falls  on  the  perimeter  of  the  net  and  the 
pole  of  this  plane  appears  at  the  origin.  With  c  vertical 
the  plane  of  the  Schmidt  net  is  the  predominant  basal 
plane  of  the  sample,  and  the  data  are  transformed  from 
the  perimeter  to  the  middle  of  the  net.  The  great  circle 
of  the  plane  of  best  fit  then  must  pass  through  the  origin. 
Accurate  visual  assessments  of  the  linear  and  planar 
preferred  orientations  of  the  sample  are  possible  on  this 
net  because  of  minimal  distortion  of  linear  distance  near 
the  origin.  With  these  planes  for  mapping,  the  impor¬ 
tance  of  net  distortion,  sample  preparation,  and  mea¬ 
surement  accuracy  are  minimized  by  providing  optimal 
views  of  the  crystal  fabric. 

Unit  vectors  in  the  Cartesian  coordinate  directions 
were  used  in  eq  1 1  to  obtain  a  unit  vector  representing 
each  crystal  in  a  sample.  These  unit  vectors  are  related 
to  an  orthogonal  coordinate  system  of  eigenvectors 
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Figure  5.  Relationship  between  the  Cartesian  and  eigenvector  coordinate  systems  for  visualizing  alternate 
Schmidt  net  representations  of  ice  fabric  data. 


(i ,  ,  i  2  ,  i  3)  by  the  matrix  P  as 

•j  =  A>ij'i  (57) 

This  relationship  is  depicted  in  Figure  5  with  eigenvec¬ 
tors  n,  c,  and  s  as  the  (primed)  unit  eigenvectors.  The 
relationship  between  a  general  vector  in  the  Cartesian 
system  and  its  transform  in  the  eigenvector  system  is 
(Hildebrand  1965) 

“m  =  />_1  “m  =  um  (58) 

and  P  is  termed  the  transformation  matrix.  Because  P  is 
orthogonal,  transformations  using  P  are  orthogonal. 
Orthogonal  transformations  maintain  length  and  pre¬ 
serve  angle,  and  can  be  interpreted  as  a  combination  of 
rotations  and  reflections.  The  order  of  the  eigenvector 
placement  in  P  determines  the  transformation,  and  ei¬ 
ther  n  or  c  is  made  vertical. 


APPLICATIONS 

We  analyzed  the  ice  fabric  of  many  samples  of  first- 
year  sea  ice  taken  from  the  Beaufort  Sea.  Both  orthogo¬ 
nal  and  dependent  variable  least-squares  analyses  were 
used  and  compared.  These  cases  are  presented  in  Figure 
6  in  an  arbitrary  order  from  highest  to  lowest  mean 
normal  error  of  the  best  plane  obtained  from  the  or¬ 


thogonal  analysis.  In  all  cases  this  approach  provided  a 
lower  bound  to  the  normal  error  of  the  dependent 
variable  methods.  The  individual  samples  are  desig¬ 
nated  by  the  dependent  variable  analysis  with  mean 
normal  error  closest  to  that  of  the  orthogonal  analysis. 
The  dashed  lines  connecting  the  mean  normal  errors  for 
the  dependent  variable  solutions  are  erratic.  Different 
dependent  variable  choices  produced  dramatically  dif¬ 
ferent  planes  with  widely  varying  normal  errors.  The 
coordinate  axis  most  nearly  normal  to  the  plane  of  the 
data  provided  the  lowest  mean  error.  Diminishing  mean 
normal  error  generally  corresponds  to  diminishing  pla¬ 
nar  spread.  However,  the  minimizations  of  these  two 
parameters  are  not  equivalent,  explaining  the  lack  of 
perfect  agreement  between  the  trends  plotted  in  Figure 
6.  The  mean  normal  errors  for  most  of  these  cases  are 
less  than  0. 1 ,  corresponding  to  planar  spreads  of  less 
than  1 2°.  These  results  indicate  that  these  samples  have 
planar  fabrics. 

The  same  cases  were  also  analyzed  for  preferred 
alignment  of  the  c-axes.  Mean  squared  normal  distance 
between  the  data  and  the  best  vector  c  is  1  -  Xy  The 
maximum  normalized  eigenvalue  and  the  linear  spread 
are  displayed  in  Figure  7.  Values  of  X'3  that  exceed  0.9 
correspond  to  linear  spreads  of  less  than  15°.  Linear 
spread  is  larger  than  planar  spread  because  it  represents 
angles  with  a  particular  line  in  the  plane,  while  planar 
spread  represents  angles  with  the  plane  itself.  Compar¬ 
ing  Figures  6  and  7  we  observe  that  increasingly  planar 


11 


Figure  6.  Mean  squared  normal  error  for  several  sea  ice  samples  from  orthogonal  and  x-t  y-  and  z- 
de pendent  variable  least-squares  analyses.  Planar  spread  from  the  best  plane  is  given  for  these  same 
cases.  The  cases  are  arranged  arbitrarily  according  to  of  the  orthogonal  analysis  and  named  in 
sequence  according  to  the  best  dependent-variable  solution. 


Figure  7.  Maximum  normalized  eigenvalue  and  linear  spread  in  degrees  for  several  sea  ice  samples.  The  order 
and  labeling  of  the  cases  is  the  same  as  in  Figure  6. 
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8, Linear  Spread  (deg.) 


Figure  8.  Schmidt  net  plots  of  cases  Zj,  zj0  and  Z/j 
including  the  intersections  with  the  hemisphere  of  the 
best  plane,  the  pole  P  of  this  plane,  and  the  best  vector 

c. 


sea  ice  fabrics  do  not  necessarily  correspond  to  the 
degree  of  preferred  alignment. 

Schmidt  net  plots  of  the  data  representing  cases  z  | ,  z  J0 
andr, ,  are  presented  in  Figure  8.  and  parameters  devel¬ 
oped  from  these  data  are  given  in  Table  1 .  Included  on 
each  net  are  the  intersections  with  the  hemisphere  of  the 
plane  of  best  fit,  the  pole  P  of  this  plane,  and  the  vector 
c  of  best  linear  fit.  As  required,  the  best  line  is  always 
contained  in  the  best  plane.  In  case  z,  the  normalized 
eigenvalues  are  approximately  equal,  the  planar  and 
linear  spreads  are  large,  and  the  matrix  A  is  nonsingular. 
These  conditions  indicate  that  the  data  are  uniformly 
distributed  over  the  surface  of  the  sphere.  Cases  z  ^  and 
z,,  are  very  different  from  z(,  but  are  similar  to  each 
other.  For  this  group  of  samples  X ,  is  smallest  for  z  M , 
indicating  that  these  data  are  the  most  closely  repre¬ 
sented  by  a  plane.  Sample  zj0  has  the  smallest  X2  value 
together  with  a  small  and  is  the  most  linear  case  in 
this  group.  The  distribution  of  data  indicated  by  the 
eigenval  ues  is  confirmed  by  the  angles  that  quantify  the 
planar  and  linear  spreads  a  and  8,  respectively,  for  each 
sample.  The  matrix  A  is  nearly  singular  in  the  latter  two 
cases  as  =  0  and  the  condition  number  is  large. 

A  group  of  samples  taken  in  close  proximity  and 
from  the  same  vertical  position  in  the  ice  sheet  should 
have  similar  structure.  Together  the  individual  and 
collective  analyses  of  these  samples  provide  quantita¬ 
tive  measures  of  comparison  that  indicate  the  sample 
size  needed  to  represent  the  ice  fabric  at  that  location. 
Three  samples  designated  zJqz  z  |B0  were  taken  from 
an  ice  core  at  1 .3  m  from  the  surface.  The  Schmidt  net 
plots  of  both  the  data  and  the  computed  fits  are  pre¬ 
sented  in  Figure  8  for  zj0  and  in  Figure  9  for  the  other 
samples  and  the  composite  data.  Each  of  these  cases  has 
nearly  the  same  planar  structure,  pole  and  linear  struc¬ 
ture.  This  similarity  is  quantified  by  the  eigenvalues, 
planar  spreads  and  linear  spreads  given  in  Table  1 . 

All  the  Schmidt  nets  presented  above  have  been 
standard,  depictingthe  lower  hemisphere  on  thexy-plane. 
Alternative  Schmidt  nets  are  obtained  for  the  hemi¬ 
spheres  below  the  best  fit  and  predominant  basal  planes 
using  eq  57  and  a  P  matrix  with  different  eigenvector 
placement  in  each  case.  As  a  result  of  these  transforma¬ 
tions  the  position  of  the  data  on  the  net  shifts,  correcting 
for  sample  preparation  error  and  allowing  a  visual 
assessment  of  the  distances  from  the  points  to  the  best 
line  and  plane.  Figure  1 0  gives  the.vy  and  the  alternative 
nets  for  sample  zJq .  The  data  are  near  the  perimeter  of 
the  net  in  the  standard  .vy-plot.  The  net  in  the  best  plane 
corrects  the  sample  preparation  error  and  displays  a 
balanced  distribution  of  data  on  opposite  sides  of  the 
net.  The  net  in  the  basal  plane  depicts  the  points  as  a 
single  group  near  the  center  of  the  net  where  linear 
distance  is  accurately  represented  and  the  linear  and 
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Figure  9.  Schmidt  net  plots  of  cases  z t0M,  zl0B  and 
composite  case  zl0  including  the  intersections  with  the 
hemisphere  of  the  best  plane,  the  pole  P  of  this  plane, 
and  the  best  vector  c. 


Figure  10.  Schmidt  net  plots  of  case  zj0  in  the  standard 
xy  -plane  with  the  z -axis  vertical,  in  the  best  plane  with 
n  vertical,  and  in  the  predominant  basal  plane  with  c 
vertical. 


Table  1.  Normalized  eigenvalues  and  condition  number  of  matrix  A,  and  planar  and  linear 
spread  of  the  data  for  selected  cases.  The  type  of  the  distribution  is  interpreted  from  these 
parameters. 


Case 

N 

h 

X 2 

Xi 

a 

(deg) 

5 

(deg) 

Cond(A) 

Type  of 
distribution 

-i 

90 

0.2940 

0.3510 

0.355 

29.3 

56.4 

1.21 

Uniform 

bi 

77 

0.00235 

0.0755 

0.922 

2.3 

13.5 

392 

Planar/linear 

rto 

44 

0.00452 

0.0417 

0.954 

3.2 

10.4 

211 

Linear/planar 

4b 

22 

0.00073 

0.0863 

0.913 

1.4 

14.2 

1250 

Planar/linear 

23 

0.00140 

0.0813 

0.917 

1.8 

14.9 

653 

Planar/linear 

2J0 

89 

0.00359 

0.0638 

0.933 

2.7 

12.5 

260 

Planar/linear 

planar  fits  can  be  readily  evaluated  and  compared  to 
other  samples. 

CONCLUSIONS 

A  quantitative  characterization  of  ice  fabrics  is  critical 
for  understanding  the  mechanical  properties  of  sea  ice, 
but  was  not  previously  available.  An  orthogonal  least- 
squares  analysis  of  uniaxial  crystal  orientation  data  was 
developed  from  geometric  arguments  with  unit  vectors 
representing  individual  crystal  orientations.  Minimi¬ 
zation  of  the  perpendicular  distances  with  a  best  line  or 
plane  provided  an  eigenvalue  problem  that  was  identical 
to  that  obtained  by  other  investigators  using  different 
methods.  Normalizedeigenvaluesgive  the  mean  squared 
normal  distance  of  the  data  from  the  line  or  plane,  and 
corresponding  eigenvectors  provide  the  dominant  c- 
axis,  planar,  and  basal  plane  orientations,  and  the  di¬ 
rection  of  columnar  crystal  elongation.  The  preferred  c- 
axis  orientation  is  always  contained  in  the  plane  of  best 
fit.  The  method  is  the  basis  of  a  relatively  simple 
algorithm  for  computer  analysis  of  large  volumes  of 
orientation  data. 

The  formulation  of  a  least-squares  method  greatly 
influences  the  results.  This  observation  was  demon¬ 
strated  for  many  samples  of  first-year  sea  ice  by  com¬ 
paring  the  mean  squared  normal  distance  of  the  data 
with  planes  obtained  using  the  classical  dependent 
variable  least-squares  approaches  and  the  orthogonal 
method.  The  dependent-variable  solutions  produced 
dramatically  different  planes  of  best  fit  with  erratic  and 
widely  varying  normal  errors .  The  error  approached  the 
minimum  given  by  the  orthogonal  method  when  the 
dependent  variable  direction  was  almost  normal  to  the 
best  plane.  Orthogonal  least-squares  and  other  analo¬ 
gous  methods  producing  the  eigenvalue  problem  all 


rely  on  the  implicit  assumption  of  equal  measurement 
error  in  all  coordinate  directions  at  all  points  on  the  unit 
sphere.  We  find  that  this  assumption  is  a  first  approxi¬ 
mation  for  optical  data  obtained  with  the  universal 
stage. 

Normalized  eigenvalues  give  the  mean  squared  nor¬ 
mal  distance  between  the  data  on  the  unit  sphere  and  the 
plane  through  the  origin  normal  to  the  corresponding 
eigenvector.  This  measure  of  the  planar  fit  to  data  is 
equivalent  to  the  variance  in  classical  dependent  vari¬ 
able  least-squares  methods.  The  majority  of  the  sea  ice 
samples  studied  had  planar  fabrics,  and  several  had 
aligned  fabrics,  characterized  by  mean  squared  normal 
distances  of  less  than  0.1  with  the  data.  These  relative 
distances  are  well  represented  by  the  angular  measures 
of  linear  and  planar  spread;  however,  the  distance  and 
angular  error  measures  are  not  equivalent.  For  these 
samples  increasingly  planar  orientations  of  the  c-axes 
do  not  correspond  to  increasingly  linear  fabrics. 

The  normalized  eigenvalues  allow  quantitative  com¬ 
parisons  between  samples,  and  of  composite  data  repre¬ 
senting  a  collection  of  several  samples.  Proximate 
samples  from  the  same  vertical  position  in  the  ice  sheet 
had  nearly  identical  structure.  This  similarity,  quanti¬ 
fied  by  the  computed  eigenvalues  and  eigenvectors, 
was  displayed  on  Schmidt  nets  forthe  individual  samples 
and  the  composite.  The  capability  to  view  the  data  and 
the  analytical  results  on  Schmidt  nets  in  the  planes 
defined  by  the  eigenvectors  was  developed,  providing 
information  to  improve  the  interpretation  of  the  data 
and  the  fits  by  minimizing  the  importance  of  sample 
preparation,  net  distortion,  and  measurement  accuracy 
limitations. 

The  eigenvalue/eigenvector  analysis  has  been  ap¬ 
plied  to  fabrics  displaying  other  patterns,  including 
multi-maxima  and  small  circle  girdles.  Fitting  these 
data  with  a  best  line  or  plane  yields  results  that  can  only 
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be  qualitatively  interpreted.  However,  when  a  fabric 
displays  several  point  maxima  it  may  be  possible  to  ob¬ 
jectively  subdivide  the  data.  The  best  line  for  each  sub¬ 
set  would  quantify  the  orientation  and  alignment  of  the 
corresponding  point  maximum.  The  intersection  of 
cones  with  orientation  diagrams  were  drawn  by  Kohnen 
and  Gow  (1979)  to  bound  small  circle  girdle  fabrics.  A 
small  circle  girdle  fabric  could  be  assessed  quantita¬ 
tively  by  finding  a  circular  cone  with  its  axis  defined  by 
the  best  line  and  its  apex  at  the  origin.  The  surface  of  the 
cone  would  be  located  to  minimize  the  sum  of  the 
squared  distances  with  the  data,  and  variability  mea¬ 
sures  could  again  be  defined  in  terms  of  distance  and 
angle. 
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